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Using Monte Carlo simulations, we demonstrate that including ionized impurity scattering in 
models of Siji.^-jGea: is vital in order to predict the correct Hall parameters. Our results show good 
agreement with the experimental data of Joelsson et.al. 
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I. INTRODUCTION 

There is much current interest in using Si/Ge strained 
layer superlattices and multiple quantum wells to fabri- 
cate optoelectronic and electronic devices - one example 
being better high frequency performance for heterojunc- 
tion bipolar transistors. In order to design effective de- 
vices we need to know the carrier concentration Nh and 
the conduction mobility /ic, which are usually obtained 
using Hall measurements. However, to extract Nh and 
/Xc from the experimental results, we need to know the 
Hall coefficient rn - but experiments in Si(i_a;)Gea; 0,0 
have shown ru has values that cannot be accounted for 
by existing attempts at modeUing (see e.g. 0,3). Also, 
the Hall factor th allows an interesting test of the accu- 
racy of the band parameters because it is quite sensitive 
to both the non-parabolicity and anisotropy of the bands. 

The reason that Hall data from existing experiments 
P, Q is not well explained by existing models is most 
probably either because (a) these models used a lin- 
early interpolated bandstructure, or (b) they omitted 
the important effects of ionised impurity scattering. Al- 
ternatives to a linearly interpolated bandstructure have 
been given by Rieger and Vogl _5j, using an empirical 
psuedopotential approach, and Hayes suggests a similar 
method based on results of a parallelized Car-Parrinello 
scheme Some other papers discussing models and 
bandstructure parameters for SiGe alloys are 0, H, 0] ■ 

Here we attack the second of these shortcomings, by 
calculating the Hall parameters using a comprehensive 
Monte Carlo simulation of hole motion in Si(i_2;)Gea; 
strained to a silicon substrate. These simulations are 
aided by the incorporation of a new optimised scheme 
for treating ionised impurity scattering di^. Our results 
show that the inclusion of ionised impurity scattering 
sucessfully explains the experimental data, leaving the 
bandstructure interpolation as a secondary effect. 

The paper is organised as follows: section Hll describes 
the Monte Carlo model used in the mobility and Hall 
calculations; section iHll compares the results of the sim- 



ulations for linearly interpolated band parameters with 
experiments. Finally in section IIVI we present our con- 
clusions. 



II. THE MONTE CARLO MODEL 

Our Monte Carlo simulations [Tlj of this system include 
a full k.p band structure calculation 12] and all impor- 
tant scattering processes: optical and accoustical defor- 
mation potential, alloy, and (in particular) ionised impu- 
rity scattering [l3l . We treat the effect of the electric 
and magnetic fields classically, giving the holes continu- 
ous trajectories in fc-space - comparison with experiment 
for p-Ge systems indicates that this approximate treat- 
ment is adequate for the field strengths we consider. We 
used a similar Monte Carlo model recently for studies 
of hot hole lasers 1^. All simulations are done at 
a lattice temperature of 300K. They follow the progress 
of a single hole through a large number of scatterings 
(typically ~ 10^ x 100), with the ergodic theorem being 
used to justify the use of the time-average as an ensem- 
ble average. Often a very large number of scatterings 
was needed in order to make sure the statistical errors 
in the Hall parameter were small enough. Our material 
parameters are listed in tabled 

The Monte Carlo simulations used the standard over- 
estimation technique where for each scattering process 
the post-scattering direction of the hole was chosen at 
random, and the differential scattering rate was overesti- 
mated by an isotropic rate just higher than its maximum 
value. Normally this would be very slow for the ionised 
impurity scattering, especially at low impurity concen- 
trations. Consequently, we use an optimised scheme 
that weights the choice of scattering angle by the angular 
dependence of the scattering rate, and thus avoid gener- 
ating a large proportion of inefficient overestimations. 



III. RESULTS 

The Hall co-efficient rn is very sensitive to both the 
shape and non-parabolicity of the bandstructure, as well 
as the anisotropy of the scattering processes. This means 
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that the role of ionised impurity scattering can be crucial 
- not only the presence or absence of impurities makes 
a difference; different impurity concentrations can pro- 
duce Hall factors that vary with alloy concentration in 
both quantatively and qualitatively different ways. As 
a result we need to compare our simulation results with 
experiment to ascertain whether our model is sufficiently 
accurate and include all the necessary important features. 

Experimental results for the Hall factor in relaxed SiGe 
have been reported by Chen et.al. for impurity con- 
centrations of about n ~ lO-'^^cm"^; and in SiGe strained 
to a Si substrate by Joelsson et.al. for n lO^^cm"'^. 
A comparison of these experimental results with Monte 
Carlo simulations has been done by Dijkstra 0)0] j whose 
simulations used linearly interpolated band parameters 
and neglected ionised impurity scattering. We have done 
equivalent simulations to those of Dijkstra, and include 
the results here. 

To calculate the Hall factor, we use the efficient esti- 
mator of the Hall mobility hh presented by 0, ^3 > 



^J■H 



V ■ X ^ 



(3.1) 



and then since rn = fin/fJ-c, the drift mobility is 
e 1 



V X B ] ■ { f X B 



(3.2) 



We did our simulations with no electric field, and a 
magnetic field in the [001] direction, in accordance with 
experiment; magnetic field strengths of i? = 0.30, 1.00, 
or 2.00T were investigated. We used a strain equivalent 
to that as if the Sii-^Ge^; alloy had been grown as an epi- 
taxial layer on a (OOl)Si substrate. This means the usual 
strain parameter Czz > 0, since the lattice parameter of 
Si is smaller than that for Ge. 

Our alloy scatteringpotential was Uq = 0.51eV. Note 
that some authors on the basis of fits to mobilities 

and other transport parameters, prefer a stronger alloy 
scattering potential (~ leV); other discussions of the ef- 
fect of alloy potential (and high impurity concentrations) 
can be found in 0, . We do not expect our results for 
th to be particularly sensitive to the strength of the alloy 
potential, as alloy scattering is isotropic. This assump- 
tion was confirmed by a set of simulations on undoped 
SiGe using the Monte Carlo program ofDijkstra||. How- 
ever, we would expect the drift mobility to diminish for 
higher concentrations of germanium is a stronger alloy 
scattering potential were chosen. 

Comparison with experimental results can be seen 
on fig. n where those for ~ 2.5 x lO^'^cm^'^ (•) and 
^ 7 X lO^^cm^^ (ir) show a Hall coefficient decreasing 
as the Ge content increases. These clearly do not agree 
with our undoped Monte Carlo simulations at 0.3T (□), 
although Monte Carlo simulations at the higher field of 
2T (A) do also have dip at low-Ge fraction. This dip at 
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FIG. 1: Hall coefficient comparison with experiment Q for 
strained Sii-^Gea; alloy, at B = 0.30T, as a function of 
Ge fraction. Monte Carlo: □ —undoped; o — 3.16 x 
10^** cm"^ ; A = undoped 2T; Experiment l]: • = 2.08- 
2.98 X 10^* cm"^; * = 6.59-7.65 x 10^** cm"^;" Error bars are 
not included because the symbols become obscured, but all 
simulation errors are listed in table HTl 




FIG. 2: Hall coefficients for strained Sii-iGe^: alloy, at 
B = 0.30T, as a function of Ge fraction. Monte Carlo: 
□ =undoped; = lO^'^cm"^; -I- = lO^cm-^*; ^= 10^** cm"^; 
o = 3.16 x lO^^cm"^; lO^^cm"^. Error bars are not in- 
cluded because the symbols become obscured, but are listed 
in table El 



2T has been previously reported for Monte Carlo simula- 
tions without ionised impurity scattering [^Q. However, 
the 2T simulations have a minimum rn at 15% Ge, with 
a subsequent increase - at 40% Ge rn is well above 1. 
This is contrary to experiments, where rn continues to 
decrease up to 36% Ge. 

However, it is clear that for undoped SiGe (fig. ^ □), 
there is no dip like that seen in experiments on doped 
material - hence, results for undoped SiGe do not ex- 
plain the experiment. Our simulations with B = 0.3T 
and Na = 3.16 x 10^*cm~"^(o) are shown on fig. ^ (also 
repeated on fig. |2Jl . These can be seen to match remark- 
ably well with the experimental data taken at Na = 2.5 
and 7 x lO^^cm"'^, making it clear that it is important 
to include the ionised impurity scattering in a model in- 
tended to reproduce the experimental value of rn- 

We can conclude that when impurity scattering domi- 
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FIG. 3: Hall coefficients for strained Sii-^^Gea; alloy, at 
B = LOOT, as a function of Ge fraction. Monte Carlo: 
□ ^undoped; = lO^'^cm"^; + = lO^'^cm"^; ^= 10^** cm"^; 
o = 3.16 X lO^^cm"^; ^= lO^^cm"^. Error bars are not in- 
cluded because the symbols become obscured, but are listed in 
tabled They are smaller for lower doping concentrations, are 
typically smaller than ±0.1, and are no greater than ±0.211. 



IV. CONCLUSIONS 

Comparison of our Monte Carlo simulation results for 
doped strained Si(i_3;)Gea; with experimental results 
have clearly shown that ionised impurity scattering has 
an important role in determining the Hall parameter rn- 
This was shown by simulations which without ionised 
impurity scattering failed to give agreement with experi- 
ment, but with it gave good agreement. This was despite 
using a linearly interpolated set of band parameters, so 
clearly the interpolation scheme is less important than 
the impurity scattering mechanism. 

Following on from this work, it would be good to re- 
duce the statistical errors in the Hall co-efficients further. 
This would enable us to make comparisons between dif- 
ferent interpolation schemes for the Si(i_2.)Ge2; parame- 
ters. Also, the behaviour of the Hall parameter in relaxed 
Si(i-a;)Ge2; would be interesting to consider. 
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FIG. 4: Hall coefficients for strained Sii_2,Gea; alloy, at 
B = 2.00T, as a function of Ge fraction. Monte Carlo: 
□ ^undoped; = lO^^cm^^; ± = lO^'^cm"^; ^= 10^*cm~^; 
o = 3.16 X 10^**cm~^; ^= lO^^cm"^. Error bars are not in- 
cluded because the symbols become obscured, but are listed in 
table^ They are smaller for lower doping concentrations, are 
typically smaller than ±0.1, and are no greater than ±0.192. 



nates, th stays low up to high Ge concentrations, a case 
which contrasts with that studied by Dykstra, who did 
not include impurity scattering. Our results thus agree 
much better with experiment than his - indeed, the im- 
portance of impurity scattering was suggested by Dykstra 
as a reason for the discrepancy he noticed. The full range 
of rn values from the Monte Carlo results are shown on 
figs. [2121 El for a range of B values. These emphasise the 
role of impurity scattering - the rn values in undoped 
and low doping simulations are significantly higher than 
those for higher doping. 
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Parameter 


Si 


Ge 


Units 


valence-band structure L 


-o.oo 


OA KO 

-30.53 


dimensionless 


M 


-3.64 


-4.64 


dimensionless 


N 


-8.75 


-33.64 


dimensionless 


spin-orbit splitting Dio 


0.044 


0.300 


eV 


deformation potential a 


2.1 


2.0 


eV 


b 


-2.2 


-2.1 


eV 


d 


-5.3 


-7.0 


eV 


do 


29.3 


40.0 


eV 


lattice constant ao 


5.43095 


5.6579 


Angstrom 


elastic stiffness cn 


16.56 


13.064 


10" dyn/cm^ 


Cl2 


6.39 


4.885 


10" dyn/cm^ 


C44 


7.95 


6.857 


10" dyn/cm^ 


mass density p 


2.328 


5.3243 


g/cm^ 


optical phonon energy Eop 


0.063 


0.037 


eV 


dielectric constant e 


11.90 


15.90 





TABLE I: Material parameters for silicon and germanium, 
as used in our simulations. They are similar to those used by 
Hinkley and Singh Hfl]. 
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Hall CO 


-efficient 


rn 


















B (T) 


Ge% 
















lo 




(cm^ 






















0.0 


15.0 


17.0 


18.0 


18.5 


19.0 


0.30 


40 


2.417 


± 


0.032 


2.056 


± 


0.052 


0.725 


± 


0.071 


0.464 


± 


0.118 


0.687 


± 


0.092 


0.761 


± 


0.161 


0.30 


30 


1.820 


± 


0.043 


1.496 


± 


0.075 


0.493 


± 


0.091 


0.358 


± 


0.132 


0.471 


± 


0.115 


0.717 


± 


0.202 


0.30 


20 


1.157 


± 


0.039 


1.017 


± 


0.101 


0.297 


± 


0.124 


0.466 


± 


0.151 


0.440 


± 


0.133 


0.494 


± 


0.196 


0.30 


10 


0.456 


± 


0.022 


0.475 


± 


0.139 


0.244 


± 


0.164 


0.240 


± 


0.189 


0.513 


± 


0.164 


0.637 


± 


0.314 


0.30 


5 


0.408 


± 


0.013 


0.253 


± 


0.160 


0.402 


± 


0.184 


0.373 


± 


0.211 


0.557 


± 


0.184 


0.715 


± 


0.349 


0.30 





0.737 


± 


0.011 


0.726 


± 


0.210 


0.601 


± 


0.196 


0.597 


± 


0.226 


0.642 


± 


0.200 


0.493 


± 


0.313 


1.00 


40 


1.788 


± 


0.016 


1.691 


± 


0.029 


0.745 


± 


0.035 


0.563 


± 


0.039 


0.531 


± 


0.076 


0.638 


± 


0.107 


1.00 


30 


1.343 


± 


0.019 


1.245 


± 


0.034 


0.551 


± 


0.044 


0.443 


± 


0.048 


0.442 


± 


0.095 


0.641 


± 


0.136 


1.00 


20 


0.903 


± 


0.022 


0.825 


± 


0.043 


0.288 


± 


0.057 


0.403 


± 


0.062 


0.413 


± 


0.120 


0.421 


± 


0.125 


1.00 


10 


0.421 


± 


0.015 


0.337 


± 


0.059 


0.289 


± 


0.077 


0.444 


± 


0.080 


0.387 


± 


0.156 


0.444 


± 


0.148 


1.00 


5 


0.423 


± 


0.011 


0.485 


± 


0.069 


0.390 


± 


0.087 


0.545 


± 


0.089 


0.509 


± 


0.174 


0.536 


± 


0.190 


1.00 





0.733 


± 


0.010 


0.767 


± 


0.074 


0.605 


± 


0.093 


0.652 


± 


0.096 


0.514 


± 


0.190 


0.716 


± 


0.211 


2.00 


40 


1.304 


± 


0.019 


1.346 


± 


0.018 


0.695 


± 


0.046 


0.514 


± 


0.052 


0.575 


± 


0.076 








2.00 


30 


0.981 


± 


0.021 


0.982 


± 


0.025 


0.596 


± 


0.064 


0.499 


± 


0.078 


0.520 


± 


0.095 








2.00 


20 


0.672 


± 


0.021 


0.661 


± 


0.029 


0.354 


± 


0.058 


0.369 


± 


0.098 


0.491 


± 


0.121 








2.00 


10 


0.374 


± 


0.016 


0.392 


± 


0.038 


0.290 


± 


0.077 


0.366 


± 


0.126 


0.517 


± 


0.157 








2.00 


5 


0.428 


± 


0.013 


0.454 


± 


0.044 


0.467 


± 


0.110 


0.461 


± 


0.141 


0.408 


± 


0.174 








2.00 





0.722 


± 


0.013 


0.779 


± 


0.046 


0.619 


± 


0.133 


0.709 


± 


0.152 


0.563 


± 


0.192 









TABLE II: Monte Carlo results for SiGe Hall Co-efficient 
vh- a full set of simulation results with data for drift and 
Hall mobilities, mean energies, time of flight, diffusion, and 
numbers of blocks required can be found in |2H| 



Drift Mobility 


B (T) 


Ge% 
















log 


10 Na (cm ^] 






















0.0 


15.0 


17.0 


18.0 


18.5 


19.0 


0.30 


40 


1517.0 


± 


3.18 


1457.0 


± 


3.26 


931.3 


± 


258.60 


478.6 


± 


0.78 


391.3 


± 


0.40 


388.0 


± 


12.14 


0.30 


30 


1156.0 


± 


3.21 


1116.0 


± 


10.23 


637.6 


± 


1.08 


390.4 


± 


0.57 


323.7 


± 


6.41 


305.2 


± 


0.52 


0.30 


20 


860.7 


± 


2.26 


827.4 


± 


1.94 


498.9 


± 


2.60 


311.2 


± 


1.52 


257.4 


± 


5.49 


238.6 


± 


0.30 


0.30 


10 


598.2 


± 


1.31 


794.7 


± 


402.90 


385.4 


± 


18.17 


240.3 


± 


0.35 


196.6 


± 


0.18 


187.9 


± 


11.53 


0.30 


5 


516.6 


± 


1.08 


502.7 


± 


1.06 


337.1 


± 


0.61 


218.6 


± 


2.44 


177.2 


± 


0.16 


162.6 


± 


0.25 


0.30 





496.2 


± 


1.20 


481.4 


± 


1.28 


324.8 


± 


0.61 


206.4 


± 


0.28 


166.1 


± 


0.16 


149.8 


± 


0.19 


1.00 


40 


1309.0 


± 


2.01 


1281.0 


± 


3.84 


791.8 


± 


2.10 


477.4 


± 


0.85 


391.3 


± 


1.09 


379.0 


± 


1.42 


1.00 


30 


1054.0 


± 


1.91 


1032.0 


± 


3.09 


633.8 


± 


1.70 


403.0 


± 


25.07 


318.3 


± 


0.89 


304.3 


± 


1.16 


1.00 


20 


816.9 


± 


1.74 


795.1 


± 


2.41 


495.6 


± 


1.36 


309.8 


± 


0.56 


286.3 


± 


64.90 


238.2 


± 


0.64 


1.00 


10 


587.7 


± 


1.19 


575.7 


± 


5.71 


376.0 


± 


1.04 


241.2 


± 


1.92 


196.2 


± 


0.55 


182.3 


± 


0.44 


1.00 


5 


510.9 


± 


1.04 


507.2 


± 


15.85 


337.0 


± 


0.96 


217.3 


± 


0.40 


177.0 


± 


0.50 


162.2 


± 


0.45 


1.00 





492.2 


± 


1.18 


478.3 


± 


1.46 


324.6 


± 


0.96 


205.8 


± 


0.39 


165.6 


± 


0.49 


149.7 


± 


0.43 


2.00 


40 


1133.0 


± 


3.04 


1111.0 


± 


3.31 


769.2 


± 


5.06 


474.3 


± 


2.23 


389.9 


± 


2.18 








2.00 


30 


946.8 


± 


2.79 


929.4 


± 


3.36 


622.3 


± 


4.73 


388.1 


± 


2.21 


317.9 


± 


1.79 








2.00 


20 


761.3 


± 


2.29 


745.5 


± 


2.71 


490.7 


± 


2.68 


311.8 


± 


6.37 


251.5 


± 


1.42 








2.00 


10 


567.6 


± 


1.70 


556.0 


± 


2.01 


374.7 


± 


2.08 


240.5 


± 


1.37 


195.9 


± 


1.11 








2.00 


5 


499.1 


± 


1.55 


488.7 


± 


1.79 


333.9 


± 


2.36 


216.7 


± 


1.26 


178.8 


± 


3.33 








2.00 





482.8 


± 


1.77 


469.0 


± 


1.75 


323.6 


± 


2.70 


205.5 


± 


1.24 


165.3 


± 


0.97 









TABLE HI: Monte Carlo results for SiGe drift mobility fic- 
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